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We study quantum interference effects due to electron motion on the Kagome lattice in a per- 
pendicular magnetic field. These effects arise from the interference between phase factors associ- 
ated with different electron closed-paths. From these we compute, analytically and numerically, the 
superconducting-normal phase boundary for Kagome superconducting wire networks and Josephson 
junction arrays. We use an analytical approach to analyze the relationship between the interference 
and the complex structure present in the phase boundary, including the origin of the overall and fine 
structure. Our results are obtained by exactly summing over one thousand billion billions (~ 10 21 ) 
closed paths, each one weighted by its corresponding phase factor representing the net flux enclosed 
by each path. We expect our computed mean-field phase diagrams to compare well with several 
proposed experiments. 

PACS numbers: 74.50. +r 



I. INTRODUCTION 

Interest in Kagome structures has been fueled by sev- 
eral recent developments. The first one refers to a cen- 
tral question in magnetism: the T = order of the 
two-dimensional (2D) nearest-neighbor-coupled Heisen- 
berg antiferromagnet. The Kagome lattice seems to 
be the first 2D spin- 1/2 model with vanishing further- 
neighbor interactions which appears to have a disordered 
ground statetl. Second, measurements of the heat capac- 
ity of 3 He absorbed on graphite at millikelvin temper- 
atures have been recently interpreted using a Kagome 
lattice structureEl. Third, measurements on the layered 
oxide SrCrg_ a; Ga4+ a ;Oi9, with Kagome-like layers, are 
currently attracting considerable attentiond. 

The connection between the ideal Kagome network and 
the real structures mentioned above (e.g., 3 He absorbed 
on graphite and SrCrg-ajGa^^Oig) is somewhat unclear 
because of the very important effects of disorder, impu- 
rities, three-dimensionality, etc., present in those materi- 
als. On the other hand, superconducting networks, made 
with electron-beam lithography, offer the possibility of 
experimentally studying, for the first time, nearly per- 
fect Kagome structures. When immersed in an extec-, 
nally applied magnetic field, superconducting nctworksu 
made of thin wires, proximity-effect junctions, and tun- 
nel junctions exhibit complex and interesting forms of 
phase diagrams (i.e., the resistive transition temperature 
as a function of the magnetic field) . 

The goals of this paper are to analytically study 
quantum interference effects of very many (~ 10 21 ) 
electron closed paths on a Kagome lattice in a trans- 
verse magnetic field B; and from these to theoreti- 
cally predict a measurable quantity, the field-dependent 
superconducting-normal phase boundary T C (B), for su- 
perconducting Kagome networks. The study of quantum 
interference effects is also important for other, physical 
processes. For instance, von Delft and Henleyu have re- 



cently investigated destructive quantum interference be- 
tween different paths that connect the same initial and fi- 
nal configurations in order to study spin tunneling events 
on a Kagome lattice. 

II. SUMMARY OF RESULTS 

The principal results presented in this paper are the 
following: (I) an analytic study of electron quantum in- 
terference effects from sums over magnetic phase factors 
on closed paths; (2) a very efficient computation of these 
"discrete path integrals" in closed-form expressions; and 
(3) from the latter, the derivation of T C {B). The lattice 
path integrals obtained here are "many-loop" generaliza- 
tions of the standard "one-loop" Aharonov-Bohm-type 
argument (where the electron wave function picks up a 
phase factor e 4 * each time it goes around a closed loop 
enclosing a net flux $). 

The calculation of lattice path integrals will enable us 
to obtain the phase boundaries^ through an iterative ap- 
proach. The spirit of the approach follows Feynman's 
programme: to derive physical quantities in terms of 
"sums over paths" . This method is considerably different 
from the standard ones for computing T C (B) (e.g., nu- 
merically diagonalizing the Ginzburg-Landau equations 
for large and/or irregular structures). We will examine 
these issues in more detail below. 



III. QUANTUM INTERFERENCE FROM SUMS 
OVER PATHS 

A. Physical interpretation 

The physics of T C {B) is determined by the electronic 
kinetic energy because the applied field-i induces a dia- 
magnetic current in the superconductor!-! This current 
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(i.e., a velocity) determines the kinetic energy of the sys- 
tem. In other words, the kinetic energy can be written 
in terms of the temperature as 



2m 



2m*£(7y 



T C (B) - T c (0), 



where, for any superconductor, m* is twice the elec- 
tron mass, and £(T) = £(0)/(l - T c ( J B)/T c (0)) 1 / 2 , is the 
temperature-dependent coherence length. We thus con- 
sider the electronic kinetic energy, 



(ij) 



exp(iAij), 



(1) 



on a discrete Kagome lattice in a magnetic field, where 
Aij = f J A-dl is the line integral of the vector potential 
along the bond from i to j. Throughout this paper, we 
set equal to one the constant factor in A^, namely 2ir 
divided by the flux quantum hc/2e. 

The lattice path integral (or moment) of order I, sum- 
marizing the contribution to the electron kinetic energy 
of all paths of Z-steps, is defined as 



mi = (Vil^lVi), 



(2) 



where denotes a localized one-site electron state. Let 
us examine more closely the physical meaning of mi . The 
Hamiltonian H is applied I times to the initial state |i/>i), 
resulting in the new state H l \ipi) located at the end of the 
path traversing I lattice bonds. Because of the presence of 
a magnetic field, a factor e tAij is acquired by an electron 
when hopping through two adjacent sites i and j. 

The above expectation value mi is non-zero only when 
the path ends at the starting site. The geometric sig- 
nificance of mi thus becomes clear: it is the sum of the 
contributions from all closed paths of I steps starting and 
ending at the same site, each with a phase factor of e 1 *^ 
where $p is the net flux enclosed by the closed path P, 
namely 



mi= ^ e 

All closed paths 



(3) 



Here, quantum interference arises because the phase fac- 
tors of different closed paths, including those from all 
kinds of distinct loops and separate contributions from 
the same loop, interfere with each other. Sometimes, the 
phases corresponding to subloops of a main path cancel. 

Several examples of different paths on a Kagome lat- 
tice and their respective contributions to the lattice path 
integrals are illustrated in Figure 1. Figure 2 shows three 
different paths traveling through the same loop and their 
separate phase factor contributions to m@ . It will be seen 
below that the phase factors describe, in a gauge invari- 
ant way, the electron interference effects, due to the pres- 
ence of a magnetic field, on the transition temperature; 
and they are the source of the rich structure present in 
the phase diagram. 



B. Analytical calculations 

We have worked out a considerable number of lattice 
path integrals (up to the 38th order) for the Kagome lat- 
tice in a uniform magnetic field. The first five moments 
can be easily computed by hand, while the higher-order 
moments can be most conveniently obtained by using re- 
cursion relations and symbolic manipulation programs. 
We have manually computed moments for the first 12 
orders. The correctness of the calculated moments is as- 
sured by the consistency of the results obtained by hand 
and by computer. The details of the computational tech- 
niques used to compute high-order moments will be pre- 
sented elsewhere. 

It is instructive to explain how the first seven lattice 
path integrals are obtained. This will also clarify their 
physical meaning. Let <j) be the flux through an elemen- 
tary triangular cell in the Kagome lattice. Since there is 
no path of one step for returning an electron to its initial 
site, mi is always equal to zero. There are four closed 
paths of two steps each [retracing each other on one bond 
(•<-»■), where the • indicates the initial site], thus 



m 2 



= 4 e 



iO<f) 



= 4 = 



where z is the coordination number of the lattice. 

There are four 3-step closed paths enclosing a triangu- 
lar cell [two counterclockwise (-v)> like the one shown in 
Fig. 1(a), and two clockwise (-v)]- Thus 



m 3 



2-v+2-V= 2e ilt > + 2e- l,p = 4cos0. 



There are 28 closed paths of four steps each: four re- 
tracing twice on one bond (•<-*■); twelve starting from a 
site connecting two adjacent bonds and retracing once on 
each bond («-►■«-►); and twelve moving two bonds away 

and then two bonds back to the original site ( — >—*)■ 
Since all of them enclose no area (i.e., no flux), then 



m 4 = 4 



+ 12- 



+ 12 



= 28. 



There are 60 closed paths (30 counterclockwise and 
30 clockwise) with five steps each. Three of these steps 
enclose a triangular cell while the other two steps retrace 

each other (e.g., -V* - *; 'V^i -^Vi •*-> Vi V"* - *'; an d 
V-^)- Thus 



m 5 



30 e^ + 30 e" 



60 cos < 



Among the 6-step closed paths, four of them enclose 
an elementary hexagonal cell (like the one shown in 
Fig. 1(b)), four paths go twice around the same trian- 
gle, 12 "hourglass" paths ( txi ) surround two triangular 
cells (6 enclosing 2<f) and 6 enclosing — 2</>), and 244 paths 
enclose no net flux. It follows then that 



m 6 



= 244 + 16 cos 20 + 4 cos 60. 
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Among the 7-step closed paths, 28 of them enclosing 
adjacent triangular and hexagonal cells (14 counterclock- 
wise and 14 clockwise) contribute 14e 7 ^ + 14e~ 7 ^ = 
28 cos 70 to 7717. Since a triangular elementary cell can 
be traversed in 756 different ways, with four steps out of 
seven enclosing no area, it follows that 

nij = 756 cos0 + 28 cos 70. 

Below we present the results for mg through 77115, while 
mi6 to m2o are listed in Appendix A and m 2 i to to 38 will 
not be presented. They are 

m 8 = 2412 + 416 cos 20 + 96 cos 60 + 80 cos 80 , 
rag = 9216 cos + 76 cos 30 + 36 cos 50 
+ 756 cos 70+ 120 cos 90, 

toio = 25804 + 7560 cos 20 + 1860 cos 60 

+ 2480 cos 80+ 100 cos 100+ 20 cos 140, 

mn = 112420 cos + 2816 cos 30 + 1276 cos 50 
+ 14608 cos 70 + 4400 cos 90 + 44 cos 110 
+ 44 cos 130 + 176 cos 150 , 

TO12 = 290956+ 1119680 cos20 + 656 cos40 

+ 33120 cos 60 + 51984 cos 80 + 4560 cos 100 
+ 36 cos 120 + 1104 cos 140 + 672 cos 160 
+ 16 cos 220, 

mi3 = 1385436 cos + 66872 cos 30 + 30680 cos 50 
+ 248872 cos 70 + 104364 cos 90 + 2756 cos 110 
+ 2184 cos 130 + 8372 cos 150 + 1456 cos 170 
+ 52 cos 210 + 156 cos 230 , 

mi4 = 3405448+ 1774248 cos20 + 31284 cos40 

+ 559160 cos 60 + 929320 cos 80+ 127876 cos 100 
+ 2632 cos 120 + 35532 cos 140 + 32256 cos 160 
+ 1960 cos 180 + 56 cos 200 + 1260 cos 220 
+ 672 cos 240 + 28 cos 300 , 

mis = 17292440 cos0+ 1307960 cos30 + 627624 cos50 
+ 3985380 cos 70 + 2047920 cos 90 + 97860 cos 110 
+ 68340 cos 130 + 246500 cos 150 + 74760 cos 170 
+ 1700 cos 190 + 3420 cos 210+ 10140 cos 230 
+ 1680 cos 250 + 120 cos 290 + 300 cos 310 . 

Note that the even (odd) moments depend only on the 
even (odd) harmonics of the flux. These moments are ob- 
tained by summing an enormous number (10 21 ) of closed 
paths, each one weighted by its corresponding phase fac- 
tor representing the net flux enclosed by the path. 

These moments, and the 77121 through 777,38 n °t shown, 
will be used in section IV to obtain superconducting- 
normal phase boundary for the Kagome network. 



C. Differences between our approach and the 
traditional moments and Lanczos methods 

In electronic structure calculations there is a method 
to compute the density of states called the moments 
method. This is similar to our approach in the sense 
that both compute moments. However, also in statistical 
analysis moments are often used. Moments are simply a 
concept of wide applicability to a large number of dissimi- 
lar problems. For instance, the half-width half-maximum 
of a probability distribution can be expressed in terms of 
its second moment; the center of mass of a distribution 
can be expressed in terms of its first moment; and the 
skewness or degree of asymmetry in terms of its third mo- 
ment. Any mathematical function (e.g., probability dis- 
tribution, Raman lineshapc, or any other function) can 
always be expressed in terms of an expansion using an 
infinite number of moments; although typically the first 
ten or so moments provide an excellent approximation 
to most functions. In our system, we have analytically 
computed moments to extremely high order (the first 38 
moments). 

We list several important differences between the stan- 
dard "moments method" and our problem. The typical 
use of the moments method: (i) focuses on the computa- 
tion of electronic density of states (instead of supercon- 
ducting T c 's); (ii) is totally numerical (instead of mostly 
analytical); (Hi) is done at zero magnetic field (instead of 
obtaining expressions with an explicit field dependence) ; 
(iv) does not focus on the explicit computation of path 
integrals; and (v) does not study the physical effects of 
quantum interference (which is at the heart of our calcu- 
lation and physical interpretation). 

In conclusion, the traditional use of the moments 
method in condensed matter is significantly different 
from the approach and problem studied here. 

Another way to diagonalize Hamiltonians is called the 
Lanczos method. This method directly obtains the tridi- 
agonal form, without computing the moments; thus dif- 
fers in a significant way from the approach used here 
(where the explicit computation of the moments is one 
of our goals, since they can be used for other electronic 
property calculations). Furthermore, it is not convenient 
to use standard Lanczos method in our particular prob- 
lem because it is extremely difficult to directly derive the 
parameters and the states of the iterative tridiagonaliza- 
tion procedure. This is so because of the presence of the 
magnetic field. On the other hand, the moments method 
provides standard procedures to diagonalize a matrix af- 
ter the moments are computed. 

Our approach has features in common with the Lanc- 
zos method and features in common with the moments 
method (e.g., the computation of the moments, which is 
absent in the canonical Lanczos technique) . More impor- 
tantly, the five differences listed above {(i) — {v)) make our 
approach to this problem quite different to the standard 
(purely numerical) implementations of both the Lanczos 
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and moments methods. 



IV. KAGOME SUPERCONDUCTING 
NETWORKS AND JOSEPHSON JUNCTION 
ARRAYS 

The complex structure in T C (B) is essentially a result of 
quantum interference effects due to the electron kinetic 
energy and the multi-connectedness of the networks in 
a magnetic field. The magnetic fluxes through the ele- 
mentary cells are useful parameters to characterize the 
interference effect. At zero magnetic field, the quantum 
interference effect is absent, and therefore the resistive 
transition temperatupe,should have a peak. 

Mean field theoryula-Q is very effective in providing 
a quantitative description of the phase diagrams. For 
wire networks, the mean field expression is given by the 
Ginzburg-Landau equation expressed in terms of the or- 
der parameters at the nodes. For a iunction array, one 
has a set of self-consistent equations!] for the thermally 
averaged pair wavefunctions of the grains. Such equa- 
tions are linearized near the transition point, and the 
highest temperature (i.e., top eigenvalue) at which a non- 
trivial solution first appears is identified as the transition 
temperature. This can be seen because the kinetic en- 
ergy of the system can be written in terms of the tem- 
perature as -h 2 /(2m*£(T) 2 ) ~ T C (B) - T c (0). There- 
fore, one is lefLln. find the top spectral edge of an eigen- 
value problemoljl3. In summary, these equations can be 
mapped into a tight-binding Schrodinger problem for an 
electron hopping on a lattice immersed in a magnetic 
field and T C (B) is determined by the kinetic energy of 
the electrons (i.e., by the Hamiltonian H). 

Now we are going to find the top spectral edge, which 
is proportional to the transition temperature. The best 
strategy for deriving eigenvalues and eigenvectors of sym- 
metric matrices is to first reduce the matrix to some 
simple form, which is typically tridiagonal. All meth- 
ods employed today (almost exclusively numerical, and 
not analytical as done here) are based on modifications 
of the Jacobi method. The latter has been in widespread 
use long time before the Lanczos and moments methods 
were developed. 

We choose a normalized initial electron state l^i), 
which is strongly localized at one site, and perform the 
following expansion: 



H\tp„) = Jn\^n-l) +Pn\lpn) + 7n+l I i>n+l ) > 



(4) 



with the condition |^o) = 0. The iJ-matrix in the basis 
\tp n ) is obviously in a real tridiagonal form. The expan- 
sion is useful because finite truncations give good approx- 
imations to the quantity we desire, i.e., the top spectral 
edge. 

Each new state in this method expands outward by 
one more step from the site where the starting state is 
located. Thus an nth-order truncation can cover a region 



of radius n on the networks. Furthermore, the param- 
eters in the truncated Hamiltonian are gauge-invariant 
quantities. 

The non-zero matrix elements, /J's and 7's, can be ex- 
actly expressed in terms of the moments of the Hamilto- 
nian by using the following novel and systematic proce- 
dure. First, define the auxiliary matrix M with the first 
row elements given by Mq^ = mi- The other rows are 
evaluated by using only one immediate predecessor row 
through 



M n _ u+2 - M„_i,i M„_i, m 



Mi 



n-1,1 



J2 M n , h M n - X 



l-k, 



(5) 



k=0 



where n and I > 1. The /3 n 's and 7„'s are obtained from 
the elements of the second and third columns as 



0n = M n . 



1.1 



and 



j n = t/M n _i 



M, 



n-1,1- 



(6) 



(7) 



Note that elements in the first column M n fi are always 
equal to 1. 

Equations (5)-(7), which are not present in Ref. 6, al- 
low a more systematic computation of the tridiagonal 
elements (3 n and 7„. Furthermore, only a few moments 
were computed in Ref. 6, and none of them applied to 
the Kagome lattice. The much higher-order moments 
computed here (involving billions of billions of paths), 
coupled with the systematic evaluation of the parame- 
ters f3 n and j n , give an unprecedented level of accuracy 
in the derivation of the structure in T c (4>). 

Below we explicitly express the first few and 7 pa- 
rameters in terms of the moments. They apply to any 
type of lattice. 



01 = 

02 = 

03 = 



0. 

m 3 

z 

m^z 2 — 2mnm^z + m 



and 



where 



7i 

72 

73 



a = UI4Z — m 2 — z 3 , 



4 



and 

A = m e m4Z — m§m\ 



m 6 z + 27715TO4TO3 + 2m 5 ni3Z 



nirz 



^>m^m\z + m^z 2 



m\ + m\. 



Using the lattice path integrals obtained, the Hamilto- 
nian matrix elements {j3 and 7) can be readily computed. 
For instance, the second order truncation of the Hamil- 
tonian is 



^/ma 



2 

2 cos< 



The presence of the cos0 matrix element corresponds 
(physically) to the Little-Parks oscillations observed in a 
single superconducting loop. Its corresponding top eigen- 
value is 




= ( COS + \J COS 2 + 16 ^ . 



The other eigenvalue is not physically important because 

(2) (2) 

it does not cross T c (</>)• T c (0) contains information on 
the quantum interference effects produced by eight paths 
(four from m 2 and four from ms). These closed paths, 

with at most three steps, cover a very small region of 

(2) 

the lattice. The corresponding AT C V (0) is the top curve 
shown in Figure 3. 

The third order truncation of the Hamiltonian is 



V 



= I 2 









y/m,2 







rru _ frn 3 \2 
mo ^ mo ' 



m± _ (223)2 _ m 



rn 2 



m 2 



ra^rn^—ra^ra^—ra 1 ^ 



2 

COS( 







7(5 



^(5- cos20)/2 
Its corresponding top eigenvalue is 
8 cos 



cos20)/2 

7 cos 0+cos 3<p 
10-2 cos 20 



Tj 3 >(0) 



where 



and 



3A 



31 - 11 cos 20 32(1 
+ — n — ~ + 



cos 2(f)) 



3A6i 



A = 5 — cos 20, 



9A 2 9* 



e = 

512 cos 3 102 cos - 22 cos 30 7 cos + cos 30 
27A 3 + 3A 2 A 

+ ^(54 cos 5 20 - 2285 cos 4 20 + 17832 cos 3 20 

- 72982 cos 2 20 + 171866 cos 20 - 179509) i 



We have analytically obtained all the eigenvalues of H^ 3 ' . 
However, the top one is the only one with physical signif- 
icance to our problem, because the two other eigenvalues 
(not shown here) do not cross it. Tc 3 \(j)) contains terms 
up to 100, and the corresponding AT C (0) is shown in 
Figure 3 (curve labeled by 3). T c (0) exhibits a few 
terms which follow the Little-Parks cos0-type oscilla- 
tions and also many additional contributions that go well 
beyond the single-loop effect since they are rigorously de- 
rived by summing 96 closed-paths on the lattice. 
The fourth-order truncation is given by 



here 



74 



2 











cos < 



a/(5-cos20)/2 



v(5-cos20)/ 2 7 -i + c :: s 2 3 / 74 

\ 74 u/v ) 



1 132 - 64 cos 20-6 cos 40 + 20 cos 60-2 cos 80 
51- 20 cos 20 + cos 40 ' 

U = 344 cos 0-109 cos 30 - 240 cos 50 
+ 570 cos 70- 124 cos 90 + 7 cos 110, 



and 



V = 692 - 449 cos 20-8 cos 40 
+ 104 cos 60 - 20 cos 80 + cos 100. 

We have also analytically obtained all the eigenvalues of 
H( 4 \ They exhibit terms which follow the cos0-type 
oscillations and also very many additional contributions 
which are rigorously derived by summing 1144 closed- 
paths on the lattice. 

In general, Tc n \(f>) is obtained from the following mo- 
ments: m 2 , m 3 , and m 2 „-i. These contain infor- 
mation on the quantum interference effects due to closed 
paths of 2n — 1 steps. 

In the following, we physically analyze how the approx- 
imations work for the phase diagrams and how the local 
geometries of the lattices affect the structures in them. 
Figure 3 shows the second through eighth, and also the 
tenth, fifteenth, and nineteenth order approximants to 
the phase diagram. The fine structures that appear in 
the lower order approximants become sharper when the 
order is increased. Eventually, the dips become cusps. 
Since the 15th to 19th approximants are essentially iden- 
tical, close convergence to the infinite system size has 
been achieved. The cusps at = 1/8, 1/4, 3/8, 5/8, 3/4 
and 7/8 are well established in the highest order shown. 
The resulting phase diagram originates from a competi- 
tion between phase coherence at different length scales, 
which can be explicitly read from the corresponding ex- 
pressions of the moments. 

Quantitatively, there are several differences between 
the T c (0) for the Kagome and, for instance, its perhaps 



5 



most similar Bravais lattice: the triangular network. The 
following simple argument illustrates a clear quantitative 
difference. For the triangular lattice, the energy is mini- 
mized when cj) — n^>Q, where n is an integer, <p is the flux 
through an elementary triangular cell, and $o is equal 
to 2ir times the quantum of flux. For the Kagome case 
this is also true: when (f> = <I>o, then <ph = 6<&o, where 
(f>h is the flux enclosed by an elemental hexagon. How- 
ever, when <fi = $0/6, the hexagons are in their lowest 
energy configuration, but the triangles are not. Thus, 
this particular value of the flux (</> = $0/6) minimizes 
the energy of the hexagonal elementary plaquettes in the 
Kagome lattice and none of the elementary cells in the 
triangular lattice. This difference will be reflected in the 
values for T c at <j) = $o/6- 

Another difference between the phase boundaries of 
the Kagome and other regular lattices (e.g., square, tri- 
angular, hexagonal) is that the latter all have a promi- 
nent minima when the flux per elementary plaquette is 
equal to $o/2, while the Kagome does not. 

This concludes our calculation and the main scope of 
our paper. In the next two sections, we briefly discuss 
experimental issues and some relations to the magnetic 
ground state of the XY Kagome Hamiltonian. 

V. HEISENBERG ANTIFERROMAGNETISM 
VERSUS SUPERCONDUCTIVITY IN KAGOME 
LATTICES 

In superconducting wire networks and arrays of 
Josephson-coupled islands, the measured T c (<p) is often 
obtained from the maximum (or some other nearby) 
value of dR/dT, where R is the resistance; i.e., the mea- 
sured T c is close, but not necessarily equal, to the real 
T c . This measured T c can be theoretically computed with 
good accuracy by solving linearized mean-field equations 
(instead of solving, for instance, the original XY Hamil- 
tonian for the Josephson-coupled islands). Mean-field is 
appropriate because fluctuations are not important as 
long as the measurements are not extremely close to the 
real T c ; which is the typical case. Also, by far the domi- 
nant effect on any given superconducting island is given 
by the nearest neighboring islands, and mean-field con- 
siders an average over them. A better approximation 
(called cluster mean-field) consists ip considering a ther- 
mal average over further neighbors! This procedure in- 
corporates a significant amount of thermal fluctuations. 
Another approach to incorporate fluctuations uses real- 
space renormalization-group techniquesH. Even the sim- 
plest mean-field for N coupled islands gives N coupled 
nonlinear complex differential equations which in general 
are difficult to solve. Linearization simplifies the mathe- 
matical treatment of the problem and is also a very phys- 
ical consideration because the density of Cooper pairs (or 
order parameter) is small near the transition. Thus, T c 
can be conveniently obtained from the ground state en- 
ergy of this linearized mean field model. 



After all these steps (i.e., mean-field and linearization), 
it is not clear that all the information present in the 
original (XY Hamiltonian) problem can be observed in 
a narrow part of phase space, i.e., close to the phase 
boundary. Even if the original information is still there, 
it is not clear to us that T c is the best way to access 
it. For example, the steps in between, from the (initial) 
full XY Hamiltonian to the (final) linearized mean-field 
equations, might reduce the symmetry of the problem; 
and some of the essential information will be lost; or still 
there but difficult to unmask in a transparent and quan- 
titative way. 

The antiferromagnetic XY model on the Kagome lat- 
tice has a degenerate ground state, and T c does not seem 
to reflect this feature because T c does not contain all the 
information that can be obtained from the full XY model 
(or the full nonlinear Ginzburg-Landau equations), but 
only part of it. Thus, it is not obvious to us how the 
degeneracy of the ground state of the antiferromagnetic 
XY model influences Tc. 

This paper does not focus on degeneracy issues, which 
are beyond the scope of this work, but on a many-loop 
version of the Aharonov-Bohm-type argument to derive 
the origin of the phase boundary features exhibited by su- 
perconducting networks. This approach involves moving 
electrons around loops threaded by a magnetic field, and 
examining the effect of this electron transport on both 
the electronic wave functions and ground state energy 
(Tc). 

VI. SUGGESTED EXPERIMENTS ON KAGOME 
SUPERCONDUCTING NETWORKS 

A measurement of AT c (<p) on a Kagome lattice would 
provide a phase boundary with features which are dif- 
ferent from the ones obtained in all the other networks 
studied so far (e.g., the absence of a prominent min- 
imum at <&o/2 and other features discussed in section 
IV) . More importantly, the Kagome phase boundary has 
many sharp dips, unlike, e.g., the relatively dull and fea- 
tureless phase diagram for the hexagonal network. It ap- 
pears to us that the dips (cusps) in the Kagome AT c (0) 
are the sharpest of all regular lattices. These sharp dips 
are predictedo to become even sharper when the mea- 
surement are made closer to the real T c . Furthermore, 
the current state of the technology would allow the con- 
struction of lattices which are far more perfect than the 
ones previously built. Thus, a direct comparison between 
measurements and our Figure 3 could be made. 

The Heisenberg and the XY Kagome antiferromag- 
nets have attracted considerable attention for the reasons 
mentioned in our opening paragraph and discussed in 
detail in references 1-3. These models exhibit geometri- 
cal frustration which prevents the establishment of long- 
range magnetic order. The classical Hamiltonians exhibit 
a complicated ground-state degeneracy, with both non- 
coplanar and planar states in the degenerate manifold. 
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In order to look for possible manifestations of the de- 
generacy and extensive entropy of the antiferromagnetic 
XY ground state, experiments could perhaps focus also 
away from T c . At lower temperatures, the imaging of 
the vortex lattice could provide a picture of the XY spin- 
state, because the phases of the order parameters would 
be arranged as the XY spin configuration of a Kagome 
antiferromagnet. At lower temperatures, the nonlinear- 
ity (e.g., quartic term in the Ginzburg- Landau free en- 
ergy) breaks the degeneracy present in the linear approx- 
imation. When lowering the temperature, the density of 
Cooper pairs increases and it is not possible to linearize 
the equations. Thus, the full Hamiltonian must be used. 

The variation of the phases (of the order parameter) 
between nearby islands will produce loop currents and 
their associated fields could, in principle, be measured 
using a variety of different probes like: (a) a scanning 
tunneling microscope, (b) a scanning atomic force mi- 
croscope, (c) a scanning Hall probe, and/or (d) induc- 
tive measurements like the ones employed by Martinoli's 
group.u Thus, information on the currents and/or flux 
configurations would be helpful to determine the phase 
arrangement (XY spins on a Kagome lattice). 
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VII. DISCUSSIONS 



The iterative approach presented above can be thought 
of as an analytical version of a hybrid between modified 
Lanczos and moments methods, both typically used in 
purely numerical procedures. It not only provides us 
with a systematic approximation scheme for the phase 
diagrams, but also a rather powerful tool for a qualita- 
tive and quantitative analysis of the structures in them 
and their physical origin. As the order of approximation 
is increased, more geometrical information of the lattice 
is included in the interference treatment, and more fine 
structures are resolved. The correspondence is such that 
we can draw a number of specific conclusions regarding 
how the dips of various sizes in T C (B) are related to the 
geometry of the underlying lattice. The gross structure 
is determined by the interference effects arising from few 
loops. The secondary dips or peaks are determined by 
the additional information provided by interference aris- 
ing from higher-order lattice path integrals (i.e., larger 
paths that explore more loops). In general, higher-order 
fine structures are due to interference among the cells of 
larger clusters (here including as many as 10 21 loops). 
The indexing of the dips (</> = n/8, where n = 1, 2, 3, 5, 6, 
and 7) in the phase diagrams can be easily read off from 
the periods of the truncated matrix elements. The sharp- 
ness of the dips is due to long range correlations, and 
cannot be accounted for by local geometries. We expect 
our predictions to compare well with future experiments 
on this system. 



7 



APPENDIX A: LATTICE PATH INTEGRALS: M 16 THROUGH M 20 

mi 6 = 41001292 + 25428480 cos20 + 918560 cos 40 + 9092352 cos 60 + 15316256 cos80 + 2863616 cos 100 
+ 109088 cos 120 + 890560 cos 140 + 956304 cos 160+ 112000 cos 180 + 5120 cos 200 + 53568 cos 220 
+ 42912 cos 240 + 2688 cos 260 + 192 cos 280 + 2752 cos 300 + 1440 cos 320 + 64 cos 380 , 

mi? = 218660460 cos0 + 23053020 cos30+ 11763116cos50 + 61614120 cos 70 + 36265760 cos90 + 2637380 cosll0 
+ 1726928 cos 130 + 5807812 cos 150 + 2339200 cos 170 + 1 13832 cos 190 + 135388 cos 210 + 385492 cos 230 
+ 1 12200 cos 250 + 2992 cos 270 + 9180 cos 290 + 23052 cos 310 + 4080 cos 330 + 340 cos 370 + 748 cos 390 , 

mis = 504967204 + 358082904 cos20 + 21570528 cos40 + 143901532 cos 60 + 240951024 cos 80 + 56411064 cos 100 
+ 3381816 cos 120 + 19345860 cos 140 + 22651848 cos 160 + 3827416 cos 180 + 256560 cos 200 + 1688400 cos 220 
+ 1610136 cos 240+ 195084 cos 260+ 17064 cos 280+ 136752 cos 300+ 105120 cos 320 + 7560 cos 340 
+ 720 cos 360 + 7668 cos 380+ 168 cos 460+ 12 cos 540, 

mi 9 = 2799941004 cos + 382148520 cos 30 + 208697444 cos 50 + 932000464 cos 70 + 603749320 cos 90 
+ 60413692 cos 110 + 38512392 cos 130 + 120342732 cos 150 + 5791 1316 cos 170 + 4436880 cos 190 
+ 4161456 cos210 + 11250660 cos 230 + 4359132 cos 250 + 246088 cos270 + 408804 cos 290 
+ 1005176 cos 310 + 303240 cos 330+ 10336 cos 350 + 30172 cos 370 + 67108 cos 390 
+ 12540 cos 410+ 1064 cos 450 + 2128 cos 470 + 76 cos 530+ 152 cos 550, 

TO20 = 6338429028 + 4999672400 cos20 + 445104760 cos40 + 2233225680 cos60 + 3683951360 cos80 

+ 1025139280 cos 100 + 87649740 cos 120 + 383654240 cos 140 + 471275520 cos 160 + 102042960 cos 180 

+ 9360020 cos 200 + 44470640 cos 220 + 46514880 cos 240 + 8109600 cos 260 + 883500 cos 280 + 4886400 cos 300 

+ 4442160 cos 320 + 595120 cos 340 + 69960 cos 360 + +438800 cos 380 + 328240 cos 400 + 26400 cos 420 

+ 2800 cos 440 + 24080 cos 460 + 12320 cos 480 + 200 cos 520 + 2080 cos 540 + 880 cos 560 + 80 cos 620 . 
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FIG. 3. Supercon- 

ducting transition temperature, ATj n '(<^>) = T c (0) — T c ((f>), 
for a Kagome lattice, versus the magnetic flux, (j>, through 
an elementary triangular cell, computed from the truncated 
Hamiltonian. From top to bottom, the orders of truncation 
are n — 2 (top curve), 3, 4, 5, 6, 7, 8, 10, 15, and 19. Note the 
development of fine structures and cusps. The convergence is 
monotonic. 



FIG. 1. Several examples of closed paths on a Kagome lat- 
tice and their respective phase factor contributions to the lat- 
tice path integrals, (a) Traversing counterclockwise a 3-step 
path on an elementary triangular cell contributes a term e 1 ^ 
to mz\ (b) a 6-step clockwise walk on an elementary hexag- 
onal loop contributes a term e -6 * 1 * to me,; and (c) a closed 
path of 26 steps with several links retraced. Its phase factor 
is e^ 16 ^ . The solid dots denote the starting (and ending) sites 
and the arrows specify the hopping directions on bonds. 



FIG. 2. An "hourglass" loop made of two triangles can be 
traversed in eight different ways, contributing three different 
phase factors to me- On the right-hand side, there are three 
possible different paths: Ti, T2 and T3. Path Ti (r2) con- 
tributes e^ + ^ = e 2 ^ (e-**-^ = e -2 ^) to m 6 . In path T 3 
the phases cancel and the contribution of this path to mg, is 
e %<t>-i4> _ ^ rp^g coe ffi c j en t s hi front of the T's are the number 
of paths with the same phase factor. 
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